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Abstract. — The general topic of the present paper is to study the conservation for 
some structural property of a given problem when discretising this problem. Precisely 
we are interested with Lagrangian or Hamiltonian structures and thus with variational 
problems attached to a least action principle. Considering a partial differential equa- 
tion (PDE) deriving from such a variational principle, a natural question is to know 
whether this structure at the continuous level is preserved at the discrete level when 
discretising the PDE. To address this question a concept of coherence is introduced. 
Both the differential equation (the PDE translating the least action principle) and 
the variational structure can be embedded at the discrete level. This provides two 
discrete embeddings for the original problem. In case these procedures finally provide 
the same discrete problem we will say that the discretisation is coherent. Our pur- 
pose is illustrated with the Poisson problem. Coherence for discrete embeddings of 
Lagrangian structures is studied for various classical discretisations (finite elements, 
finite differences and finite volumes) . Hamiltonian structures are shown to provide co- 
herence between a discrete Hamiltonian structure and the discretisation of the mixed 
formulation of the PDE, both for mixed finite elements and mimetic finite differences 
methods. 
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Introduction 

Many problems arising in various fields (such as physics, mechanics, fluid mechan- 
ics or flnance) are described using partial differential equations (PDEs). Although 
explicit solutions are not available in general, important classes of PDEs do present 
strong structural properties: classical examples are symmetry properties, maximum 
principle or conservation properties. It is quite essential for the numerical methods 
to provide a translation of these structural properties from the continuous level to 
the discrete level so enforcing the numerical solutions to obey qualitative behaviours 
in agreement with the underlying physic of the problem. 

Two fundamental notions arising in classical mechanics are Lagrangian and Hamil- 
tonian structures. Lagrangian systems are made of one functional, called the La- 
grangian functional, and a variational principle called the least action principle. 
From the least action principle is derived a second order differential equation called 
the Euler-Lagrange equation, see e.g. [1]. The Lagrangian structure is much more 
fundamental than its associated Euler-Lagrange equation: it contains informations 
that the Euler-Lagrange equation does not. An important example is the change of 
coordinates. The Lagrangian structure is independent from change of coordinates, 
whereas the associated Euler-Lagrange equation may completely change of nature 
(from linear to non linear for instance). A range of numerical methods forget about 
the Lagrangian to focus on the Euler-Lagrange equation itself. This is for example 
the case for the flnite difference methods. 

In this paper we address the first following question: consider a PDF deriving 
from a Lagrangian and a least action principle. When discretising this PDF, how 
is embedded the attached Lagrangian structure at the discrete level ? More precisely 
we ask whether the discretised PDE can be seen as deriving from a discrete least 
action principle associated with a discrete Lagrangian structure. Basically, in case 
the Lagrangian structure is embedded at the discrete level, then the variational 
property of the original equation (at the continuous level) may be preserved by the 
discrete problem. 

Aside from the Lagrangian systems are the Hamiltonian systems: under suitable 
conditions a Lagrangian system can be reformulated as a Hamiltonian one. PDEs 
deriving from an Hamiltonian structure inherit strong properties such as energy 
conservation or existence of first integrals. Such PDEs moreover concern a wide 
class of problems: like mixed formulations and saddle point problems [3] that usually 
are associated with an important class of numerical methods (mixed finite elements 
e.g.). 



The second question we address in this paper is whether Hamiltonian structures 
can be embedded at the discrete leveh do numerical solutions of PDEs deriving 
from a Hamiltonian structure also derive from a discrete Hamiltonian structure ? 
Again, in case a Hamiltonian structure is embedded at the discrete level, strong and 
interesting properties might be inherited by the numerical solutions. 

The concept of coherence is introduced to answer these two questions. Let us 
consider a problem having a Lagrangian structure. That is to say consider a La- 
grangian functional £ on a functional space and the attached calculus of variations. 
The solutions to the continuous problem satisfy a least action principle that reads 
a second order PDE: the Euler-Lagrange equation. Discretisation can be performed 
at two different levels: 

— either discretise the Euler-Lagrange equation by defining discrete analogues to 
the differential operators in this PDE. We will call discrete differential embed- 
ding this procedure, 

— or discretise the Lagrangian structure by defining a discrete Lagrangian func- 
tional Ch- This second procedure is called discrete variational embedding. A 
calculus of variations then can be developed at the discrete level on Ch to define 
a discrete least action principle. This process of obtaining a discrete counter- 
part of the Euler-Lagrange equation with the help of the discretisation of the 
Lagrangian functional and the discrete least action principle is also called vari- 
ational integrator. There is a wide range of work by Lubich [17] , the group of 
Marsden J18j on this topic preserving structures for ODEs. 

In case the two discrete differential and discrete variational embeddings do define 
equivalent discrete problems we will say that we have coherence. This is enunciated 
in saying that the following diagram commutates: 

discrete Lagrangian Ch 

disc. L.A.P. 

Euler-Lagrange equation '■ '■ > discrete Euler-Lagrange equation 

where L.A.P stands for least-action principle. The same notion of coherence can be 
defined relatively to Hamiltonian structures. 

In case of coherence, the discretisation firstly preserves the variational structure 
of the problem so inheriting interesting properties (such as independence with the 
coordinate system). It secondly may also preserves algebraic properties from the 
differential operators within the Euler-Lagrange equation. 

Based on this notion of coherence, the present work is an attempt to inter- 
pret numerical methods as variational integrators for PDEs deriving from a La- 
grangian/Hamiltonian structure. We will focus on a canonical example of such a 
problem: the Poisson equation. This classical problem being well documented both 
at the continuous and at the discrete levels, it provides an appropriate test case 
to improve the understanding of discrete embeddings for Lagrangian/Hamiltonian 
structures and of the embedding of variational properties from the continuous to 
the discrete level. Let us point out that such structures are not available for every 
common problems, such as convection problems or the Navier-Stokes equations. The 
structural numerical difficulties (numerical diffusion, stabilisation e.g.) attached to 
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these problems have to be hnked to this lack of variational structure. A general alter- 
native to build new numerical schemes (or variational integrators) for such problems 
would be: 

1. construct a non classical variational structure for the considered problem as in 

il^IID], 

2. embed this structure at the discrete level. 

Although this paper is mainly concerned with improving the understanding of dis- 
crete embeddings of variational structures on a classical example, it is aimed to help 
towards the development of new numerical schemes for it helps in controlling the 
second step above. 

The outline of the paper is as follows. Part |T] gives an informal introduction to 
embedding formalism. We introduce the notions of discrete differential and discrete 
variational embeddings of a given problem, as well as the concept of coherence 
between these two embeddings. 

Part HI] deals with Lagrangian systems. In section H] are defined the Lagrangian 
structure and the associated calculus of variations for fields. The Lagrangian struc- 
ture for the Poisson equation is recalled in section |5l In section [6] is presented the 
discrete embedding of a Lagrangian structure: the principle of coherence is defined 
and conforming finite element methods are here shown to satisfy the coherence prin- 
ciple for general Euler-Lagrange PDEs. 

In part mil we focus on the Poisson equation and on the coherence of two classical 
numerical methods for this problem: finite differences and finite volumes in sections 
[7] and E] respectively. 

Part II VI is concerned with Hamiltonian structures and mixed formulations. Hamil- 
tonian structure and the associated calculus of variations are presented in section [9l 
One recover the mixed form of the Poisson problem through its Hamiltonian formu- 
lation. The discrete embedding of Hamiltonian structure and the notion of coherent 
embedding are presented in section [101 coherence is shown to be naturally fulfilled 
by conforming mixed finite element methods. Coherence of mimetic finite difference 
methods, |3], for the mixed Poisson problem is analysed in the last section. 

The general notations are detailed and listed at the end of the paper. 



PART I 
DISCRETE EMBEDDINGS 



The formalism of embeddings has been initiated in [8] and further developed in 
[6], [9], |llj . [2]. This part introduces, on an informal way, the notion of discrete 
embeddings in the case of partial differential equations. The definitions will be 
detailed in partsUT t Illll and HVl For a complete introduction to embedding formalisms 
for PDEs we refer to [7]. 



1. Discrete differential embeddings 

Let fi C M"' be a bounded domain and let M denotes some functional space on Vt. 
In general, a PDE can be written as: find n G M so that, 

(1) P{u) = 0, 

where P is a differential operator. This differential operator can be written in 
different forms. As an example, the classical Laplacian operator 

1=1 

can be rewritten as 

A = divoV. 

We will adopt the following informal definition of a discrete differential embedding 
of PDE (d]): 

Definition 1 (Discrete differential embedding). — A discrete differential 
embedding of equation (QP is formally obtained by: firstly providing a discrete 
space Mfi (of finite dimension) in which the discrete solution Uh will be sought and 
secondly by providing a discrete operator Ph on M^ analogue to P. 
The discrete differential embedding of equation (J\) becomes: find Uh G M^ so that 

(2) PhM = 0. 

As an example, let us formally consider Ah, divh and V^ discrete analogues of 
the operators A, div and V (precisions being given in the sequel). This leads to 
two discrete differential embeddings for the Poisson problem: find Uh € Mh so that 
either, 

^hUh = 0, 

or, 

div h{VhUh) = 0. 

These two discrete problems do not coincide in general. Indeed, recovering the 
algebraic properties of the original differential operators (here A = div o V) at the 
discrete level (here A^ = div/j oVh) is a full problem by itself. 

A second problem related to discrete differential embeddings is that it is based on 
the form of the operator P, which form depends on the considered coordinate system 
(used to write the equation). 
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2. Discrete variational embeddings 

We here assume that equation ([T]) derives from a variational principle. Precisely 
that it is obtained from a least action principle on a functional 

namely, that the solutions of equation ([T]) coincide with the critical points of C 
Being considered a subset V G M referred to as the set of variations and denoting 
DC{u){v) the Gateau derivative of £ along a direction v eV ,we precisely have: for 
any u G M: 

P(m) = ^ DC{u){v) = {),\/ V eV. 

This is referred to as a variational formulation for problem ([T]) and the correspon- 
dence between the solutions of the PDE and the critical points of the functional a 
variational principle. 

Rather than discretising the differential operator P (i.e. to perform a discrete 
differential embedding as previously defined) we here are led to the following alter- 
native: discretise the functional C and define the discrete solutions Uh through a 
variational principle on £/,. 

Definition 2 (Discrete variational embedding). — Assume that the PDE (QP 
has a variational formulation given by a functional C and a set of variations V C M. 
A discrete variational embedding of equation (Cp is formally obtained by: firstly 
providing a discrete space Mh (of finite dimension) and a discrete set of variations 
Vh C Mh and secondly by defining a discrete functional Ch '■ Mh -^ M analogue of 
C. 

The variational embedding of equation (J\) is then defined as: find Uh G Mh such 
that 

(3) DCh{uh){vh) = Q, ^Vh^Vh, 

(i.e. Uh is a critical point of Ch for the variation set Vh). 

Discrete variational embeddings essentially possess interesting properties. The 
main one being that the underlying object supporting the discretisation does not 
depend on the coordinate system, on the contrary of the differential form ([T]) of the 
problem and therefore of its discrete differential embedding. Discrete variational 
embedding corresponds to the so-called variational integrators studied for ODEs by 
Lubich [17j and the group of Marsden 



3. Coherence of differential and variational embeddings 

As formally introduced here, there is no reason for the two discrete differential and 
variational embeddings of a given problem to provide equivalent discrete problems. 
This question is addressed considering the concept of coherence introduced in [S]. 
An informal statement is the following : 

Definition 3 (Coherence). — Let problem ([T]) have a variational formulation. 
Two given discrete differential and variational embeddings for ([T]) are said to be in 
coherence in case the two discrete problems ([2]) and ([3]) are equivalent. 



In case of coherence, the discretisation of problem ([T]) not only preserves its vari- 
ational structure but also the algebraic one of P. 

A general raised question then is: Can we find conditions ensuring the coherence 
between the discrete differential and variational embeddings ? 

In the next three parts we study the coherence for discrete embeddings of problems 
having a Lagrangian or Hamiltonian variational formulation. It turns out that one 
cannot set apart the coherence from the algebraic properties of P^ inherited from 
the one of P: more precisely for properties of integration by part type. A deeper 
insight into this relationship is gained by considering the Poisson problem. When 
performing a discrete differential embedding for the Poisson problem involving both 
a discrete gradient and a discrete divergence (to define P^), coherence is obtained in 
case these two discrete operators fulfil some discrete analogue of the Green- Gauss 
formula. This is the case for finite differences with formula ( ITTl) as detailed in remark 
m This also is the case for finite volumes with formula fl22|) in remark |5l Considering 
the mimetic finite difference (see section [TT]in part lIVp . such a discrete Green-Gauss 
formula is a priori imposed in order to build the scheme and coherence here again 
is satisfied. 



PART II 
DISCRETE VARIATIONAL EMBEDDING OF LAGRANGIAN 

SYSTEMS 



In the first section, we recall classical results about Lagrangian calculus of varia- 
tions. Section |5] deals with the Lagrangian formulation of the Poisson problem. The 
notions of discrete variational embedding and coherence in the Lagrangian case are 
specified in section [Ul As a first example, we discuss the case of conforming finite 
element methods. 



4. Lagrangian calculus of variations 

In this section we define basics on calculus of variations applied to vector fields. 
For more details, we refer to the books of Evans |12] . Giaquinta and Hildebrand 
Let us note M = L^{Q). 



Definition 4- — An admissible Lagrangian function L is a continuous function 

L-.nxRxW^ — y R 

{x,y,v) h-^ L{x,y,v) 

such that L is of class C^ with respect to y and v. Let us consider a subspace Dom(£) 
of H^{fi) satisfying: for all u G Dom(£), L(-,m(-), Vm(-)) G L^{fl). The Lagrangian 
function L defines the Lagrangian functional C: 

C : Dom(£) -^ M, 

u I — 7- / L{x,u{x),Vu{x)) dx. 
Jn 
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We are interested in vanishing the first variations of the Lagrangian functional L 
on a space of variations V . As in |14] . we could give general notion of extremals 
and variations. We take the following definitions of the notions of differentiable 
functional and extremal for C 

Definition 5 (Differentiability). — We consider a space of variations V C 

Dom(£). The functional C is differentiable at point u G Dom(£) if and only if the 

limit 

C(u + ev) — C(u) 
lim 

exists in any direction v E V. We then define the differential DC{u) of C at point 

u as 

^^/ N/ N , C(u + €v) — C(u) 
V eV ^ DC{u){v) ■= lim ^ '- ^. 

e-S>0 e 

With the above definition of differentiability, one recovers the usual definition of 
the differential in case V = Dom(£) and DC{u) is continuous on Dom(£). The 
definition given here suffices to introduce extremals: 

Definition 6 (Extremals). — A function u G Dom(£) is an extremal for the 
functional C relatively to the space of variations V C Dom(£) if C is differentiable 
at point u and: 

DC{u){v) = for any v ^V . 

Proposition 1. — The Lagrangian functional L is differentiable at point u G 

dL dL 

Dom(£) if: x i— ?■ — — (a;, n(x), Vn(x)) and x i— )■ ---{x,u{x),'Vu{x)) are in M and 
oy ov 

]\^d respectively. In such a case the differential is given for any v G Dom(£) by: 



(4) DCiu)iv) 



— — (x, u(x), Vm(x)) V + — — (x, m(x), Vu(x))-Vv 
oy ov 



dx. 



Proof. — Using a Taylor expansion of L at the point (x, u + ew, V x{u + ev)) in the 
variables y and v leads to: 

dL dL 

L(x,u + ev,V(u + ev))= L(x,u,Vu) + ev —-(x,u, Vu) + 'Vx(£v)-—-(x,u,Vu) + o(e) 
^ oy ov 

Integrating over the domain Q gives: 

Ciu+ev) = C{u)+ I ev —-(x,u(x),\/u(x))dx+ / (e'Vv)-—-(x,u(x),\/u(x))dx+o(e) 
Jn oy Jn ov 

leading to (i!). D 

Extremals of the functional C can be characterised by an order 2 PDE, called the 
Euler- Lagrange equation given in following theorem. 

Theorem 1 (Least action principle). — Let us assume the Lagrangian func- 
tional C is differentiable at point u G Dom(£) and that u is an extremal for a 
given space of variations V . Assume moreover that ^(-j u(-), Vm(-)) G H^^^iVi) and 
that the subspace Vq = {v & V,v = Q on dVt] is dense in M. Then u satisfies the 
generalised Euler- Lagrange equation: 

dL ( dL \ 

(5) T^i^i u(x), Vuix)) — div I t:— (x, u(x), Vuix)) ) = 0. 

dy \ dv J 



Proof. — Following (J4]) and using the Green formula gives: V f G Vq, 

p r f) r / f) r \ ~ 

/ —-(x,u(x),\/u(x)) — div{—-(x,u(x),\/u(x))] vdx = 0. 
Jn [dy \dv J 

which implies i^ by density of Vq in M. D 

5. Lagrangian structure for the Poisson problem 

We consider the following Poisson problem on Q for a homogeneous Dirichlet 
boundary condition: find a solution u to 

, . j -Am = / in fi, 

^^ \ M = on dn. 

For a data / G M, the problem ([6]) has a unique solution m G i/^ fl -ff^Kfi). 
For a function / G M, we consider the (admissible) Lagrangian function L: 

L{x,y,v) = -v-v- f{x)y. 

The associated Lagrangian functional C is defined on Dom(£) = H^(r2) by: 

C{u) = / -\'Vu{x)\'^ — f{x)u{x)dx, 
Jn 2 

and is differentiable on H^{Q). 

Theorem 2. — The solution of the Poisson problem ^ is the extremal u eV for 
C with the space of variations V = Hq{Q). 

Remark 1 (Generalisation). — Consider the general elliptic problem: 

- div(aVn) = / in Q, 



(7) 



u = on dQ. 



for the tensor x G ^2 i— > a(x) G R''^'^ (strongly elliptic, positive, symmetric, and 
bounded). The solution of problem ^ is the extremal u G Hq{Q) of the Lagrangian 
functional C associated with 

L{x, y, v) = -{a{x)v) ■ v - f{x)y, 

and for the space of variations V = Hq{Q). 



6. Coherence for discrete embeddings of Lagrangian structure 

In this section, we suppose given a conformal mesh Th of the domain Q. A precise 
definition for a mesh (which is not needed here) is given in definition [TTl To the 
functional spaces M and V are associate discrete counterparts Mh and Vh. We 
suppose given a discrete Lagrangian functional Ch defined on Mh. 

Remark 2. — A general definition for a discrete Lagrangian is not obvious and 
greatly differs from a numerical scheme to another. In the following, we show that 
a discrete Lagrangian for the general Euler-Lagrange equation ([5]) is available in the 
case of conforming finite element methods. Part l//il deals with the definition of a 
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discrete Lagrangian for the Poisson problem in the case of finite volume and finite 
difference methods. 

Definition 7 (Discrete least action principle). — A discrete function u^ G 
Mh is an extremal of Ch relatively to a space of variations Vh if 

(8) DCh{uh){vh) = for any Vh E Vh- 

Equation ((HI) is called the discrete Euler-Lagrange equation and the procedure 
discrete least action principle. 

Definition 8 (Coherence principle). — A discretisation procedure is said co- 
herent if the discrete differential embedding of the Euler-Lagrange equation ([5]) is 
the discrete Euler-Lagrange equation ([8]) coming from the discrete least action prin- 
ciple. This is equivalent to the commutation of the following diagram: 



L.A.P. 



disc. L.A.P. 



u solution of PDE '■ '■ > Uh solution of PDEf^ 

(E.L. equation ^) (discrete E.L. equation ^) 

In other words, the procedure is coherent if the direct discretisation of the Euler- 
Lagrange equation given by the discrete differential embedding leads to the same 
solutions as the one given by the associated discrete variational embedding. 

Coherence for conforming finite element methods. — We consider here the 
Euler-Lagrange PDE ^ together with a homogeneous Dirichlet boundary condition 
M = on dVt. For simplicity we assume that we can set Dom(£) = HQ(r2). We 
consider M^ a classical conforming finite element discretisation of H^ {VL) relatively 
to the mesh Th- For instance, in case the mesh is a triangulation, Mh denotes the 
P^'-Lagrange finite element space of continuous functions polynomial of order k in 
each cell of the triangulation (see e.g.[5l, I16j ). We then define Vh as the subset of 
Mh made of the functions Uh € Mh that vanishes on dVt. 

The discrete differential embedding of PDE @ given by the finite element method 
is: find Uh G Vh such that, for all Vh € Vh, 

(9) / (^-—{x,Uh,Vuh)vh + -^{x,Uh,Vuh) ■ Vvhjdx = 0. 

A discrete variational embedding is provided by defining the discrete Lagrangian 
Ch as: 

Ch := £|Mh- 
This definition holds since for conforming finite element methods Mh C Dom(£). 

Theorem 3 (Coherence). — The conforming finite element methods for PDE ([5]) 
is coherent. Precisely the singular points of the discrete Lagrangian Ch with respect 
to the space of variations Vh and the solutions of the discretisation ^ coincide. 

Remark 3. — By definition conforming finite element methods are based on some 
weak formulation. This weak formulation might or might not be a variational for- 
mulation of some least action principle: this is typically the case for the Poisson 
equation but not for the for convection- diffusion or Navier-Stokes equation. If the 
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weak formulation coincides with a variational formulation, the conforming finite el- 
ement discretisation is naturally coherent with the underlying variational principle 
and the associated Lagrangian structure. 



PART III 
COHERENCE OF CLASSICAL DISCRETE EMBEDDINGS 



In section [6] of the previous part we showed a first example of coherent discrete 
embedding of Lagrangian structure. In this precise case, several facilities where 
available: the discrete solution had a clear definition as a function Uh '■ ^2 — ^ M so 
that differentiation and integration had the same sense at the discrete and at the 
continuous levels. As a result the definition of a discrete Lagrangian C^ was obvious 
and natural: Ch was the restriction of C to some functional space of finite dimension. 

Such facilities are not always available: they rather are restricted to conforming 
finite element methods. In general, differentiation and integration have to be re- 
defined at the discrete level to provide a definition of a discrete Lagrangian. In this 
part we give two examples of discrete embeddings for a Lagrangian structure: finite 
differences and classical finite volumes. Coherence is proved for the discretisation of 
the general Euler-Lagrange PDE ([5]) in the case of finite differences and in the case 
of the (isotropic) Poisson problem for finite volumes. 

7. Finite differences 

We consider in this section the general Euler-Lagrange PDE ([5]) together with a 
homogeneous Dirichlet boundary condition m = on dQ. We refer to [19] concerning 
the finite difference method. For simplicity, the domain will here be set to ^2 = [0, l]"^. 
We consider a Cartesian grid 7^ of ^2 with uniform size h = 1/N, A^ G N*, in each 
direction. The results of this section can be extended to more general domains and 
more general lattices without difficulty. 

We will use the following notations here. For j = {ji,...,jd) € N"', we de- 
note j^* (resp. j^* ) the shift of j by +1 (resp —1) in the z*'* component: j^* := 
(ji, . . . , ji_i, ji+1, ji+i, . . . ,jd), J denotes the subset of N^ J := {j G N"', < j < N}. 
The point of coordinates {jih, . . . ,jdh) is denoted Xj G M'^. The considered grid Th 



then is the grid with vertices {xj, j G J}. 

7.1. Discrete differential embedding. — Let us consider the two spaces of 
mappings: S := {uh '■ Z*^ — > M} and V := {Fh '■ 'L^ — > M*^}, whose elements are 
respectively denoted (mj) and (Fj). We introduce the discrete gradient and discrete 
divergence operators: Vh : S — )■ V and div/^ : V — > S as: 



^ij>i -U\\ d r 



(10) (V,«,)j = M : , (div,F,)j = ^ 



Mj>d — Uj 



h 



where [X]j stands for the i*'* component of the vector X G M.'^. 

We consider the discrete functional spaces M^, and Vh defined as: M^ := {uh G 
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S, Mj = if j ^ J} and Vh := {uh G M^, Mj = if j G dJ}. The resort to S and V 
is motivated by the simphfications it provides in the notations only. 

With these definitions, we have the following discrete Green-Gauss formula: 

(11) V M, e V^, V F, G V : ^Fj ■ {VhUh)j h^ = -^^(div.FOj u^ h". 

jGJ jeJ 

Remark 4- — Up to now, different choices have been made for the definition of 
discrete differentiation operators. The discrete gradient is defined using a forward 
Euler scheme whereas the discrete divergence uses a backward Euler differentiation 
formula. Other choices are possible: for instance forward on the divergence and 
backward on the gradient, or either a centered scheme for both operators. The cru- 
cial point is that whatever are their definitions, the two discrete operators have to 
commute in a discrete Green-Gauss formula of type ( fTTl) to eventually ensure the 
coherence. 

Definition 9. — The discrete differential embedding given by the finite difference 
discretisation of the Euler- Lagrange PDE ([5]) reads: find Uh G Vh such that, 

(12) ( — (x, Uh, VftU/,) - div/i ( — (x,u/„ Vh.u/,) j j =0, 

for all i such that Xj G f2 (i.e. < j < N) and with: 

f9L^ ^ A dL 

(^-(x,n„V.n,)J_ = - 

dL . ^ A dL . ,^ , , 

— (X, Uh, VhUh) 1 = ^(Xj, Mj, [VhUh)i) 

7.2. Discrete variational embedding, coherence. — 

Definition 10 (Discrete Lagrangian). — We consider the following definition 
for the discrete Lagrangian functional Ch '■ Mh — > M, associated to the continuous 
Lagrangian functional C : M — ;■ M given in definition [^ 

(13) £h{uh) = Y. ^(^J' ^J' i^hUh)i) h''. 

The discrete Euler- Lagrange equation then is: 

(14) Find Uh G Vh such that : DCh{uh){yh) = for any Vh G Vh- 

Theorem 4 (Coherence). — The finite difference discretisation of the Euler- 
Lagrange PDE i^ is coherent in the sense with definition\^ Precisely: the solution 
of the discrete Euler- Lagrange equation JHM and the solution of the discretised 
Euler-Lagrange equation flT^ coincide. 

Proof. — Let us consider a solution to (TT^ . We then have: 
jeJ ^ jeJ 



( -^{x, Uh, VhUh) ] = ^rX^v ^j' {^hUh)-i> 



Using the discrete Green-Gauss Formula (fTTj) we exactly recover the discretised 
Euler-Lagrange equation ( TT2i) . D 
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8. Finite Volumes 

We consider in this section the classical finite volume method (as presented e.g. 
in |13] ) for the Poisson problem (I6j). The following definitions and notations will be 
needed. 

Definition 11 (Mesh Th, cells K, faces £ and vertices M) 

A cell (or control volume), generically denoted K, is a polygonal/polyhedral open 
subset K G Q. 

A mesh Th of the domain Q is a collection of cells K partitioning Vt in the following 
sense: 

UKeThK = n , Ki,K2eTh => either KiH K2 = (!) or Ki = K2 . 

A face (or an edge) e of some K eTh such that e C dVt is called a boundary face. 

The set of boundary faces is denoted Sq. It satisfies: dQ = Uee^-oC For every e G Sq, 

there exists a unique K gTh satisfying e C A' fl dVt: one writes e = K\dVL. 

The internal faces set £i associated with Th is the set of all geometrical subsets 

e = Ki n K2, Ki,K2 G Th and Ki 7^ K2, having non-zero {d — 1) ~ dimensional 

measure. For every e G Si, there exist a unique couple Ki,K2 G Th satisfying 

e = KiH K2: one writes e = Ki\K2. 

The faces set associated with Th is given as S := SoUSi. It provides a partitioning of 

Ux^Th^K, in the same meaning as earlier: Ue6£:e = Uk^Th^K and the overlapping 

of two distinct faces either is empty or of zero {d — 2) — dimensional measure. 

Eventually, the set of vertices associated with Th is denoted M : it contains exactly 
all the vertices of all the cells K G Th. 

Let e G S such that e C dK for K G Th. We denote n^.e the unit normal to e 
pointing outward of K. 

One shall also denote by \0\ the measure of a geometrical object O according to 
its dimension: taking d = 3, \K\ is the volume of the cell K, \e\ the area of e E S 
and \xy\ the length between two points x and y. 

Two sets of points are introduced: cells centres {xK)KeTh ^^^ boundary faces 
centres {xe)ee£o- They satisfy: 

(15) W K eTh, We e So : Xk e /i, Xe G e . 

\/e e Si : e = Ki\K2, [xK^.XKi] -L e, 

(16) Ve G £"0 : e = K\dVL, [x^, Xk] ± e. 

Conditions OlSp -f TTBj) are referred to as admissibility conditions. They impose a 
strong constraint on the mesh Th- Non conformal meshes for instance cannot fulfil 
such a constraint. 
Distances {de)e<^£ across the faces are defined as follows: 

Ve = Ki\K2 G Si : df, = IxkiXkzI > 
Ve = K\dQ E Sq : de = \xKXe\ ■ 

8.1. Discrete differential embedding. — To the mesh Th is associated the 
discrete functional space Mh- 

(17) Mh := {uh G M, Uh=^ UkXk, Uk G M} , 
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where xk is the characteristic function of the cell K: Mh is the space of piecewise 
constant functions on the cells K E Th- An element Uh G Mh will be simply denoted 
by its component: Uh = {uK)KeTh- 

The following interpolation operator I is considered: 

(18) I: feM^f eMh. fj^:=±. f fdx, W K e %. 

To Uh G Mh, are associated its numerical fluxes along the edges e E S, using a finite 
difference scheme as follows: 

(19) We = K,\K,e£. : F,,k, = !f^^i_^f^, 

y e = K\dQeSo : F,K = -^- 

de 

Definition f lTUj) of the fluxes on the boundary is the discretisation of the homogeneous 
Dirichlet boundary condition in (jU]) considered here. Of particular importance is the 
following continuity relation: 

(20) We = K,\K2e£i, F,^K, = -Fe,K, ■ 
This property motivates the following definition: 

W e E S , Fe := |-Fe,x| , for any K E Th such that e C dK . 

Definition 12. — A discrete Laplace operator Ah : Mh — )■ Mh is defined as fol- 
lows. For Uh G Mh: 

\/K eTh , {AhUh)K ■= TT7\ y^ Fe,K\e\ . 

A ^ — ' 

The discrete differential embedding given by the finite volume scheme for problem 
([6]) reads: 

(21) find Uh E Mh so that : —Ah.Uh = f^ ■ 

Remark 5. — It is interesting to consider the discrete Laplacian A/i as the compo- 
sition of the flux operator in (IT9l) with a discrete divergence operator (of finite volume 
type). Being given a flux distribution F = {Fe^K)ee£, eeax, the discrete divergence is 
given on each cell K eTh by: 

divxF = — ^ Fe,j^|e|. 

' ' e&dK 

In this way Ah in definition [TH expresses a flux balance on each cell K mimicking 
in a discrete context the divergence formula. 

Moreover we have the discrete Green-Gauss formula for a flux distribution F 
satisfying the continuity property ( l20l) and for Uh E Mh: 

(22) 

^ {([iykF)uk\K\ = - ^ Fe,K^ — ''—, ^|4||e|+ ^ Fe,i^Mi^|e|. 

K&Th e=Ki\K2e£i ^ e=K\dn££o 



15 



8.2. Discrete variational embedding, coherence. — In the continuous case, 
the energy term | Vwp is part of the Lagrangian functional C In the framework of 
finite volume method, no proper discrete gradient is available. Actually the numeri- 
cal fluxes definition ( 1T9|) involves a discrete differentiation in the normal direction to 
the mesh faces. Thus only the normal component of some discrete gradient on the 
mesh faces is approximated in the finite volume framework and not the tangential 
one. 

The discrete differentiation along the face normal direction is used to define a dis- 
crete Lagrangian functional: 

Definition 13. — The discrete Lagrangian functional Ch '■ M^ — )■ R zs defined as: 
(23) Mun e M, , £,K) = ^ J^i^f \e\de - Yl ^k^k \K\ . 

Theorem 5 (Coherence). — The finite volume discretisation for the Poisson 
problem (E]) is coherent in the sense of definitions^ Precisely: let us consider the 
discrete least action principle associated with ( 123|) .■ find Uh G M^ such that, 

Then this problem has a unique solution which is moreover the unique solution of 

m- 

Note that here Mh = Vh, the boundary condition being encoded directly in the defi- 
nition of the discrete Laplace operator (through the fluxes definition (fT9ll ) and not 
in the choice of the variation space. 

Proof. — With definition ( l23l) . Ch is clearly strictly convex and has a unique mini- 
mum. Differentiating ( l23l) . it comes: V Uh,Vh € Mh, 

DCh{uh){vh) = - ^ fj^VK\K\ + ^ --^VK\e\ 



KeTh e=K\dneeo 



e 



+ 2^ 2 K^K^-VK.,)\e\. 

e=Ki\K2€£i " 

Using the discrete Green-Gauss formula (1221) with the flux distribution F equal to 
the flux of Uh given by ( fT9l) we simply get: 

DCh{Uh){Vh) = - J2 {fK + dwKF)vK\K\ =-J2 {fK + {^hUh)K)vK\K\. 
KdTh KeTh 

Thus, the condition DCh{uh){yh) = for all Vh G Mh exactly reads —AhUh = f^. □ 



PART IV 

HAMILTONIAN CALCULUS OF VARIATIONS AND MIXED 

FORMULATIONS 



In this part let L : fi x M x M"^ be an admissible Lagrangian function as defined 
in section |H We recall the link between Hamiltonian and Lagrangian systems. 
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Moreover we try to understand the mixed formulation as a Hamiltonian system. As 
an application, we consider the Poisson problem. 
Let us recall that M = L^{n) and X = [L'^{n)f. 



9. Hamiltonian and mixed formulation 
9.1. Hamiltonian formulation. — 

Definition 14 (Legendre property). — We say that L satisfies the Legendre 

dL 

property if the mapping v H- -;:—(x, y, v) is a bijection on M.'^ for any x ^ Q, y ^W. 

ov 

If L satisfies the Legendre property, the following function g : ^2 x M x R'^ — )■ R'^ 
is well defined: 

v = g{x,y,P) ^tth p= —{x,y,v). 

Let us consider p := t^(x, y, v) as a new variable, then, 

ov 

p= —{x,y,g{x,y,p)) and g{x,y, —{x,y,v))= v . 

Definition 15 (Hamiltonian). — Let L satisfy the Legendre property. The 
Hamiltonian iJ:r2xRxR°'— j-R associated to L reads: 

H{x,y,p) =p- g{x,y,p) - L{x,y, g{x,y,p)). 

We introduce two different definitions for the Hamiltonian functional % : 
Dom(?^) — )■ R associated to H with domain Dom('H) C M x X.- 

1. Primal Hamiltonian. 

(24) 'H{u,F):= f F-Vu-H{x,u,F)dx. 

Jn 

2. Dual Hamiltonian. 



(25) 



n{u,F) := / -div{F)u- H{x,u,F)dx. 
Jn 



For both the primal and the dual cases, we formally define Dom('H) as a subspace 
of M X X on which the Hamiltonian functional is well defined. For the primal case 
we have Dom('H) C H^(r2) x X whereas for the dual case Dom(7/) C M x fldiv(^)- 

Proposition 2. — We consider a space of variations VxW G Dom('H) . Following 
the definitions in section [^ the Hamiltonian functional "H is differentiable at point 
(u, F) G Dom(H) if 

PjTJ f)]-f 

-T—(x,u,F)eM and -T—(x,u,F)eX. 
ay op 

In such a case we have: 

— In the primal case: 
(26) 

Vu- ^—(x,u,F)] +Vv-F -v^—(x,u,F) dx. 
op J oy 



D'H{u,F)-{v,G) 
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In the dual case: 



D'Hiu,F)-{v,G) 



BJ-f f rtJ-f 

div(G)M - G ■ ^— (x, M, F) - t; ( div F + ^—{x, u, F 
op \ oy ' 



dx. 



Definition 16 (Extremals). — Let us consider a space of variation V ^ W d 
00111(7/). We say that (ti, F) G Doni('H) is an extremal for l-i relatively to V x W 
ifH is differentiable at point (u, F) and: 

Dn{F,u)-{G,v) = y {v,G)eVxW. 

Theorem 6 (Hamilton's least action principle). — Let (u, F) G Dom('H) be 
an extremal for Ti relatively to V x W . Assume moreover that: 

— in the primal case: F G H^i^^Q), Vq = {v ^ V,v = on dQ} is dense in M 
and W is dense in X , 

— in the dual case: u G H^{VL), V is dense in M and Wq = {G G VF, G ■ n = 
on dQ} is dense in X . 

Then {u, F) is a solution of the Hamiltonian system: 

f)J~f 

divF = — -—{x,u,F) 

Vn = —{x,u,F). 

Proof. — Let us consider the case of the primal definition of the Hamiltonian func- 
tional "H. Since F G H^i^{Q), using the Green formula in ( 126|) gives: V(f , G) G 
VxW, 

-(divF + ^—(x,u,F))v + G(Vu- ^—(x,u,F))]dx + / vF-nds = 0. 
in dy ^ ^ dp '} Jg^ 

The boundary integral vanishes for t> G Vq. We recover (127|) by density of Vq in M 
and of W inX. D 

Corollary 1 (Lagrangian and Hamiltonian formulations) 

The solutions {u, F) of the Hamiltonian system ( 1271) are exactly the solutions of 
the Euler- Lagrange equation ([5]) under the condition 

dL 

F:=^(a;,«,V,n)). 

9.2. Application to the Poisson problem. — We consider problem (|7]). We 
recall that the Lagrangian function associated with this problem is 

L{x,y,v) = -{a{x)v) -v- f{x)y. 

dL 
The Legendre property is clearly satisfied by L since -;—(x,y,v) = a(x)v. We 

ov 
introduce the new variable p = a{x)v and the function g is given by g{x,y,p) = 

a~^{x)p. A Hamiltonian for the Poisson equation is then given by 

(28) H{x,y,p) =p- {a~^{x)p)-L{x,y,g{x,y,p)) = -a~^{x)p-p + f{x)y. 
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Proposition 3. — The Hamiltonian system ( |271) associated with ( !28l) is the mixed 
formulation of the Poisson problem, (JTj); 

- div F = / 
(29) { 



The mixed problem (1291) together with the homogeneous Dirichlet boundary con- 
dition u = on dVt has the classical mixed weak formulations: 

1. Primal formulation: look for (m, F) G HI{VL) x X such that 

\ («-^F - Vn, G)o,f7 = , VGgX 

2. Dual formulation: look for (ti, F) G Af x iJdiv(^) such that 

r (divF + /,^)o,n = 0, Vt;GM 

\ («-'F,G)o,n + KdivG)o,n = 0, VGG/Jdiv(fi) 

Proposition 4- — A solution for flHUl) zs an extremal {u,F) G -f^o(f^) x X /or i/ie 
primal definition ([21]) o/T/ relatively to the space of variations V xW = Hq{Q) x X . 
A solution for ( l3T|) zs an extremal (m, F) G M x ilfdiv(f^) /or t/ie (ifia/ definition 
( l25l) o/?/ relatively to the space of variations V x W = M x i7div(^)- 

10. Discrete Hamiltonian and Coherence 



Let us consider a discretised version of the mixed problem ( 1271) obtained by some 
numerical scheme. If the discrete problem is also an extremal of a discrete Hamil- 
tonian functional the considered numerical scheme is said to be coherent. It is the 
analogue of definition [8] to the Hamiltonian functional. Let us precise here the 
notion of coherence regarding Hamiltonian structures. We assume we have two dis- 
crete spaces Mh and Xh and an equation on Mh x X^ that we refer to as the discrete 
differential embedding for PDF flTTj) . On the other hand assume that we have a 
discrete Hamiltonian Tih '■ M^ x X^ — )■ M. If the solutions {uh,F^) of the discrete 
differential embedding for system of PDEs (1271) are the solutions of the discrete least 
action principle DTihiuhjF^) = then the discretisation is said coherent. This can 
be summarised by saying that the following diagram commutes: 

(«, F) G M X X h^ n{u, F) "''"■ """ "'"'^-> (n,„ F J E M^ x X,, ^ Hh{uh, F 



L.A.P. 



h, 

disc. L.A.P. 



(m, F) solution of PDE ([27]) '^''^- ^'^- '"'^•) (m/,, FJ solution of PDE/, 

(Hamiltonian system) (discrete Hamiltonian system) 

Remark 6 (Coherence for conforming mixed finite element) 

We point out here that the mixed finite element schemes naturally are coherent. 
Consider some classical discretisation Mh and Xh of M and X (for example Mh = 
P'^il'h) andXh = RT^ilh) the Raviart- Thomas elements of order 0, see e.g. |16j ). A 
discrete differential embedding of (127]) given by the conforming mixed finite element 



" 


dH , 


f OH 


-diY{GJuh- 


-G, ■ g^{x,Uh,Fj- 


- Vh div F^ + ^ 
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method relatively to a homogeneous Dirichlet boundary condition reads: 
look for {uh, F^) G Mh x Xh such that for any (u/j, G^) G Mh x Xh, 

(x, Ufe,Fj) dx = 0. 
In 

A discrete variational embedding is provided by defining the discrete Hamiltonian 

Tih '■ Mh X Xfi -^M. as the restriction ofH to M^ x X^. This definition makes sense 

when considering conforming finite elements where Xh C iJdiv(^)- 

So defined, the discrete embeddings naturally are coherent. 

Remark 7. — Let us note that the definition of the Legendre property is not clear 
on the discrete level. As an example, considering again the Raviart- Thomas elements 
of order for the mixed Poisson problem, we do not get Yh = Vuh, but some more 
complicated formula depending also on the source term f. 



11. Mimetic Finite Differences 

We consider the mixed Poisson problem (!29|) together with a homogeneous Dirich- 
let condition m = on dfl. The dual weak formulation (IHT]) is adopted. The spaces 
of variations are set to V^ = M = L'^{Q), and W = X = H^i^^Q). The ambient space 
M X X is equipped with the following scalar products: 

(m, v)j^j = uv dx , (F, G)j^ = / a^^F ■ G dx. 
Jn Jn 

A flux operator Q : H^{Q) C M — > X is introduced: Qu = aVu. This operator is 
adjoint to the divergence operator div : W G X — > M in the following sense: 

yueH^in), VFgI^: (F,^n)^ = -(divF,n)j^^ . 

In the Mimetic Finite Differences (MFD) framework, two discrete functional spaces 
Mh and Wh are introduced allowing to define discrete divergence divh '■ Wh — >■ Mh- 
Scalars products on Mh and Wh are then introduced. A discrete flux operator is 
eventually introduces as (minus) the adjoint of the discrete divergence. We refer to 
[3j for the MFD discretisation of diffusion problems. 

11.1. Discrete differential embedding. — A mesh Th of the domain Q is con- 
sidered, the definitions and notations in definition [TT] are adopted here. The discrete 
space Mh is set to be the space of piecewise constant functions on the mesh cells 
K & Th defined in (IT71) . We consider the notations introduced in section [8] and the 
interpolation operator I : u E M — y u^ G Mh defined in ( ITSl) . We consider on Mh 
the scalar product induced by the Euclidean structure on M: 

(32) \/uh,Vh G Mh : K,f/i]jv/^ = ^ ukVk\K\. 

KeTh 

Fluxes are associated to each face e E £. To define such fluxes, a crossing direction 
has to be considered on the faces: to every cells K E Th and to every e E £ such that 
e C dK we define n^ e the unit normal to e pointing outward from K. Numerical 
fluxes Fx,e £ K are associated to each face e in the direction n^^e- We introduce the 
continuity condition: 

(33) We = K,\K2 E Si ■■ F^^e + ^x^.e = 0. 
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The discrete space Wh is defined as 

Wh := {{FK,e) ioT K eTh, eeS s.t. e C dK that satisfies dM])} . 

The dimension of Wh is equal to the number ^S of faces in S. However, no canon- 
ical isomorphism between Wh and M*^ is available: to define such an isomorphism 
necessitates to prescribe a crossing sense on each face e E S. 
The following interpolation operator is considered: 

FeW^F' eWh : F^, := ^ /"f ■ iiK^eds. 

\e\ Je 

Definition 17. — The discrete divergence divh : Wh — > Mh is defined as, 
^FheWh, ^KeTh, (div, Fh)K := 4^ E ^^M- 

' ' ee£,eCdK 

For all F G W , we have (div F) = divh (F^) , meaning that the following diagram 
commutes: 

W -^ M 



Wh -^^ Mh 

The definition of a scalar product on Wh is not obvious. Let us consider K E Th 
and denote W^ the restriction of Wh to K. We suppose that a cell scalar product 
[■,-]k is given on each W^ G Th and thus defines the scalar product on Wh as: 

(34) yFh, GhEWh-. [Fh, GhW ■■= Y. t^'^' ^'^l^' 

the definition of [■, ■]k will be detailed later on. Relatively to the scalar products 
([321) ^'^d (IMj) . the discrete fiux operator Qh : Mh — > Wh is defined as (minus) the 
adjoint of the discrete divergence Qh = — divj^, where -k stands for the adjoint. 

Problem ( 12^ is then discretised as follows, [1]: find Uh G Mh and Fh G Wh such 
that, 

-AiWhFh =f' 

Fh - QhUh = 
which is equivalent with: find Uh G Mh and Fh G Wh such that, 

, . f y VhE Mh, [dwhFh,Vh]M^ = -[f^,Vh]M^ 

l^GheWh, [Fh, Ghlw, + [uh, div, Ghhi, = 



A way to define the elemental scalar product (1341) is to introduce a lifting operator 
TZk : W^ — y [L\K)f and then to define: 

[Fh,Gh]K := / a-\x)nK{Fh) ■ nK{Gh)dx. 
Jk 

Consistency is ensured by imposing the three following conditions on the lifting 
operators. For all i^ G 7, and all Fh G W^^: 

Wees, eCdK : nK{Fh) ■ rvK,e = FK,e ] div (7^^(F0)) = div^ (F^ ; 
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and for all constant vector c G M'^: 

TZk {c') = c. 

11.2. Discrete variational embedding, coherence. — We consider the dual 
definition (l25l) for the Hamilton functional % on M ^W relatively to H defined in 
(EHD. It reads: 

U{u, F)= f a"^F • aVu -\ f a'^F ■ F - /" f{x)u 
Jn 2 Jq Jq 

= (F,^«)^-i(F,F)^-K/)^. 

Definition 18 (Discrete Hamiltonian) . — We therefore consider the following 
definition of the discrete Hamiltonian functional: T-Lh '■ Mh x Wh — > K; 

(36) ?^,(«,, F,) := [F,, ^,w,]^^ - ^ [F,, F,]^^ - [m,, /^] ,^^^ . 

Theorem 7 (Coherence). — The MFD discretisation for the mixed Poisson prob- 
lem ( l29l) zs coherent in the sense of sectionlWi Precisely: The discrete Hamiltonian 
Tih in ( l36l) /ias a unique singular point {uh, Fh) G Mh x PF/i such that 

DHhiuh, Fh){vh, Gh) = VK, G;,) G M;, X W^. 
This singular point is the unique solution to flH3|) . 
Proof. — Differentiating Tih gives: 

D'Hhiuh,Fh){vh,Gh) = [Fh,GhVh]w^ + [QhUh,Gh]^^^ - [F/„G/,](^^ - [/^%]^^ ■ 
A singular point {uh, Fh) for Tih then satisfies: 

\/ Vh^ Mh, [Fh,GhVh]wh = [f^,Vh]M^, 

^GhEWh, [Fh,Gh]w,-[GhUh,Gh]w, =0 

which exactly reads (!H^ by using ^/, = — div^. 

One can easily prove that (1351) has a unique solution. D 

Notations 

— Q open bounded domain fi C M'^, 

— dQ domain boundary, 

— M scalar function space on Q, M = L^(f2), 

— X vector function space on ^2, X = L'^{QY, 

— L Lagrangian function L -.Vl xW^ xW^ — > M, 

— C Lagrangian functional C : M — y M, 

— H Hamihonian function H -.nxR"^ xR'^ — ^ M, 
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— "H Hamiltonian functional Ti : M x X — ;• M, 

— h index for all discrete objects. 

— Th mesh of Q, 

— £ mesh faces (or edges in dimension 2), 

— £q boundary faces, 

— Si internal faces, 

— M vertices set, 

— K one cell (or control volume) of Th, 

— e one face (or edge) of Th, 

— I "I object measure (according to its dimension), 

— H"^{Q) denotes the Sobolev space of order m, 






1/2 



for any / G H^{Q), its gradient is given by V/ 



9x2 f 



\ d.J J 
we refer to V^/ to emphasise that the derivative is respect to the x parameter, 

d 
the divergence of a vector F is div F = N^ d^^Fi, 

i=l 

HaU^) = {v e {L\n)Y, divt; G L\Q)}. 
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